
AIAA 96-1677 

A Study of Wave Propagation in a Duct 
and Mode Radiation 

F. Farassat 

NASA Langley Research Center 
Hampton, VA 

M. K. Myers 

The George Washington University 
Hampton, VA 


2nd AIAA/CEAS Aeroacoustics Conference 

May 6-8,1996/State College, PA 


For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics 
370 L'Enfant Promenade, S.W., Washington, D.C. 20024 




A STUDY OF WAVE PROPAGATION IN A DUCT 
AND MODE RADIATION 


F. Farassat* 

Fluid Mechanics and Acoustics Division 
NASA Langley Research Center 
Hampton, Virginia 

M. K. Myers** 

The George Washington University 
Joint Institute for Advancement of Flight Sciences 
NASA Langley Research Center 
Hampton, Virginia 


Abstract 

In this paper we discuss two problems of classical 
duct acoustics: (1) wave propagation in infinite ducts 
with uniform flow based on a graphical approach, and 
(2) detection of mode radiation from a duct by an 
external circular microphone array. In ( 1 ) we show that 
the wave number vectors for a given flow Mach number 
form an ellipse whose center and shape depend on the 
Mach number only. We construct graphically the 
upstream and downstream wave number vectors of 
modes that propagate in the duct. We then show how 
one can infer from this graphical approach many known 
results in duct propagation such as the mode cut-off 
concept, the direction of energy propagation, the angle 
of the radiation lobe at peak directivity using Rice’s cut- 
off ratio concept, and other qualitative results. In (2) 
we give the mathematics behind the experimental 
detection of mode radiation from a duct by a circular 
array of microphones whose axis coincides with the 
engine axis. Since the external microphone array does 
not introduce additional noise sources inside the engine, 
as does a rotating rake of microphones positioned at the 
inlet, this measurement technique may be preferable to 
use of a rotating rake. The simplicity afforded by the 
lack of sophisticated rotating parts is an additional 
advantage over the rotating array. 

Introduction 

Acoustic wave propagation through a flow in an 
infinite duct and sound radiation from a duct inlet are 
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old problems of acoustics. Both are relevant to the 
design and understanding of the noise characteristics of 
ducted fan engines. Noise is generated inside the engine 
by various flow related mechanisms, e.g., rotor wake- 
stator interaction, and then propagates through the inlet 
and exhaust ducts. Although the duct length in all 
engines is finite in length, the study of acoustic wave 
propagation in infinite ducts of circular or annular 
shapes has been very useful to engine designers both 
qualitatively and quantitatively. The basic paper on 
infinite duct acoustic wave propagation is by Tyler and 
Sofrin 1 . There are many useful results in this paper and 
elsewhere on infinite duct propagation that one should 
understand in order to explain some aspects of engine 
noise. A real engine, however, has nonuniform flow 
within the duct which complicates the problem of 
acoustic wave propagation considerably 2 3 . One would 
then like to know what modes are generated inside and 
propagated out of the inlet and the exhaust. One way of 
finding modes that escape out of the duct is by using a 
rotating microphone array 4 . However, the rotating 
microphone array produces a wake that may interact 
with a fan rotor and produce additional noise. It would 
be preferable if one could use a noninvasive 
experimental method of detecting escaping modes in 
place of a rotating microphone array. 

In this paper we intend to present two results 
related to the acoustic problems discussed in the above 
paragraph. First, we give a graphical representation of 
wave propagation in an infinite duct carrying a 
uniformly moving fluid. We show that many known 
results can be derived easily by this graphical method. 
The method relies on construction of an ellipse in wave 
number space whose shape is dependent on flow speed. 
The propagating waves as well as upstream and 
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downstream wave number vectors are then constructed 
graphically. The direction of energy propagation, the 
mode cut-off concept and the approximate radiation 
angle for each mode can all be constructed graphically 
also. Second, we propose a new method for 
experimental detection of modes escaping from a finite 
duct using an external circular array of microphones 
whose axis coincides with the engine axis. We give the 
theory behind mode detection here. Such an array has 
been constructed at Langley Research Center and will 
soon be tested in an anechoic room using a model 
ducted fan engine. 


p ' = A m J m [kj (m, n) r] exp i [ w t - m 0 - k^m, n) x] ( 4 ) 

where the k f (m,n) are obtained as the roots of 

jl[k I (m,n)R] =0 < 5 > 


numbered in consecutive order by n=l, 2, 3, . . .. Here 
J m (*) is the Bessel function of the first kind of order m. 
Let p 2 =l-M 2 . Then k a (m,n) is given from Eqs. (3) and 
(4) by 


Graphical Approach to Wave Propagation 
in a Uniform Duct With Flow 


k,(m,n) = ± [ -M ± yi-[pk t (m,n)/k] 2 


Consider a duct of uniform circular or annular 
cross section carrying a uniform flow at Mach number 
M<1. The differential equation and the boundary 
condition for propagation of small pressure perturbations 
for the hard wall case is 


■AI-m./'-wL] <,) 

in which m=0, 1,2,.... We have defined the cut-off 
ratio p mn as 


— — +M—j p'-V 2 p' = 0 
c 0 dt dx j 

(l-a) 

p - — !L^ 
Pk,(m,n) 

(7) 

=0 (on the wall) 
dr 

d-b) 

It then follows that the mode (m,n) 
p mn > 1 and decaying if P ran < 1 

is propagating if 


in which x is the axial coordinate along the duct. If 
cylindrical polar coordinates (r, 0, x) are used, and if it 
is assumed that a complex solution for the acoustic 
pressure p' exists in the form 

p / =P(r,e,x)e‘"* (2) 

then the complex amplitude P satisfies the equations 


In preparation for the graphical approach, we write 
Eq. (6) as follows 


(k, +M/P 2 ) 2 
1/P 4 



( 8 ) 


where k 1 = k a /k and k^l^/k. For use later we note 
that Eq. (8) can also be expressed as 


(1-M 2 ) — + ^P-2ikM^+k 2 P=0 (3 _ a) 

dx 2 dx 

— =0 (on the wall) (3_b) 

3r 

in which k=oi/c 0 and is the two-dimensional 
Laplacian in polar coordinates. 

Let k,(m,n) and k r (m,n) be the axial and radial 
wave numbers for the mth circumferential and nth radial 
modes. Considering a circular duct of radius R now, 
the solution of Eq. (1) for the mode (m,n) is 


+ k| = ( 1 -Mk,) 2 

Equation (8) is an ellipse in the variables (kpi^) with 
center at (-M/p 2 , 0) and with semi-major and semi- 
minor axes 1/p 2 and 1/p, respectively. The ellipse 

depends only on M. It intersects the kj axis at 
kj = 1/(1 +M)>0 and kj = - 1/(1 -M)<0. It always 

intersects the kj axis at kj = ±l. Figure 1 shows the 
ellipse described by Eq. (8) for flow Mach number 
M=0.8. Note that for M=0, Eq. (8) gives a circle of 
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unit radius in the kjk^ plane with its center at the 
origin. We now use this ellipse to study some aspects 
of wave propagation in a circular duct. 


u' = - 


V 


k,p' 


ip 0 c 0 k(l Mkj) p q Cq( 1 Mk,) 


-5, +a 2 


( 11 ) 


it Propagating and Decaying Modes : Figure 2 
shows how one determines the propagating and decaying 
modes. On different vertical axes, each corresponding 
to a circumferential mode m, plot the solutions of Eq. 
(5), i.e. k r (m,n)/k for n=l,2,. . .. Draw horizontal lines 
as shown in this figure for m=0. If these lines intersect 
the ellipse, the mode (m,n) is a propagating mode. 
Otherwise the mode is decaying. In general, we get two 
values for k a which we denote as k a+ and k a _. In Fig. 2 
we have shown k a+ (0,3)/k and k a .(0,3)/k. Note that k a 
determines the axial phase speed of the mode (m,n). 
The k a+ mode is generally called the downstream mode, 
but any propagating k a4 mode for which k/k>l actually 
moves in the negative x direction; this is the case with 
k a ,(0,3) on Fig. (2). The relationship of the phase speed 
of the mode to the velocity of energy propagation in the 
mode will be discussed below. 

jj) Mode Cut-off and Energy Flow : When the cut- 
off ratio P m „ is equal to unity, k^l/p. Since the 

semi-minor axis of the ellipse in the k,^- plane is 1/p, 
the axial wave number corresponding to 

k^l/P is kj =k a /k = -M/p 2 . Figure 3 shows the 

wave number vector for the cut-off condition at M=0.8. 
It indicates that the k a+ wave propagates upstream. To 
discover what is special about the wave number vector 
at cut-off we consider the acoustic energy flux vector 
for the modal solution given by Eq. (4). 


in which e, is the unit vector along x and u^ is the 

projection of u / on the cross-sectional plane of the duct. 
Thus, for the single-mode complex solution given by 
Eq. (4) we have 


W =Re 


p' 

1-Mk, 



P 2 k, +M 

Po c 0 ( 1 -Mk 1 ) 


/_ -»/ 
P «l +u 2 


J 


( 12 ) 


Here Re[ ] denotes the real part of [•]. We see from Eq. 
(12) that at the cut-off condition, k : = - M/p 2 , the axial 

component of the acoustic energy flux vanishes: W is 
normal to the wall and no energy propagates along the 
duct axis. 


This cut-off of axial energy flux can also be 
exhibited graphically on the wave number ellipse, 
although only approximately for a circular duct (the 
interpretation is exact for ducts of rectangular cross- 
section). For this purpose we use the large argument 
approximation of the Bessel function 




N nk r r 


cos(kjr - tl» m ) 


(13) 


with <|i ni = 7i/4+mn/2 to approximate p' by 


For isentropic disturbances in a uniform flow at 

velocity V the acoustic energy flux vector can be 
written as' 


W=(p / +p 0 V-u / )(u / + - 
=(p' + p 0 c 0 M-u')(a' 



Po c o 


+ -£^M) 
Po c o 


(10) 


In Eq. (10), u' is the acoustic particle velocity 
corresponding to p', p 0 is the density of the undisturbed 

flow, and M is the flow Mach number vector V/c 0 . 
The (primed) disturbance quantities are the actual (real) 
acoustic pressure and velocity. For an acoustic field 
represented in the complex form of Eq. (4) the 
linearized equation of linear momentum gives 


r 

2 \ 


2 

/ exp i [cot -m0 -k x -k r r - t|r n 


expilwt-me-k.x + ^r-MjiJ} (14) 


We see that, for sufficiently large k,r at any fixed 0, p' 
can be considered to be the sum of two progressing 
waves in (x,r) whose phases are 

4> ± = Mt-k 1 x*k I r = k(c 0 t-k 1 xTk 2 r) (* 5 ) 

The unit vectors normal to these planar phase surfaces 
are n a = -Vi|> t /|V<l>J, or 

_ (kpikj) (kj, ±1^) 

n = — = (lo) 

^ 1 “Mk. 
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where we have used Eq. (9). Now if we invert Eq. (16) 

to express the wave number vector k 1 s (k l ,±k 2 ) in 
terms of n ± , we obtain 

=• - ( 17 ) 

k t = fi t /(l + M -n t ) 

If we use Eq. (17) in Eq. (15) we find 

<t> =u)t- kP,r ( 18 ) 

* 1 +M n i 

where r is the position vector (x,r) in the plane of 
constant 0. The phases in Eq. (18) are precisely those 
of a pure plane wave propagating in the directions h ± in 
the (x,r) plane. The phase speeds of these waves are, 

from Eq. (18), c o ( 1 + M n ± ) . 

If the approximation is viewed locally and 

variations in the geometric factor r' 1/2 are ignored, then 
Eq. (14) implies that at points sufficiently close to the 
duct wall at least the higher order n modes can be 
interpreted as being a sum of two interfering plane 
waves propagating in the directions n ± . For each of 

these separately the relation p o c 0 holds so that 

the energy flux associated with each is given from Eq. 
(10) by 

„ _/2 

W = — — (1 +M -h)(c 0 h + V) (19) 

Po c o 

The last factor in Eq. (19) is the energy propagation 
velocity, or the group velocity V Q , of the plane wave: 

V Q = 0 ^ + V =c 0 (n +M) 

and it is this quantity associated with the approximate 
plane wave representation of the duct modes that is 
convenient to consider in conjunction with the wave 
number ellipse. 

If we again refer to Fig. (3), we see that at cut-off 
the component of V G in the x direction is 

V G *e t =c 0 ( -cos0 c +M) =c 0 ( -M + M) =0 (21) 


Here we have used the fact that, as seen from Fig. 3, 
tan 0 c =p/M and thus cos 6 C =M, where 0 C is the angle 

that k at the cut-off condition makes with the kj -axis. 
Equation (21) means that the acoustic energy flux vector 
is normal to the wall and no energy is propagated along 
the axis. 

To find the direction of group velocity in our 
graphical construction follow the vector diagram of Fig. 

4. This diagram is in the kjkj plane so that M must be 

scaled up to M/p 2 to get the direction of the group 
velocity as shown in this plane. From the figure, it is 
clear that at the cut-off condition the group velocity 
points upward, and thus the energy flux vector is normal 
to the wall. 

iii) The Angle Of The Main Lobe Of Radiation 
From An Inlet : Rice, Heidmann and Sofrin 6 have 
proposed a rule for finding the angle of the main lobe 
of radiation from an inlet where the duct Mach number 
is M and the inflow Mach number (flight speed) is M b . 
By our graphical method we can give a simple 
interpretation of their rule as follows. Intuitively, it is 
obvious that the phase velocity of a radiating mode at 
the inlet should be based on M, the duct flow Mach 
number. But the angle of the main lobe of radiation is 

in the direction of the group velocity using M m instead 

of M in Fig. 4. Figure 5 shows this construction. Here 
0' is the sought angle. The mathematical reasoning is 

as follows. Remembering that kj<0, we have, from 
Fig. 5: 

tan0' = ^ ^ — — (22) 

k, + M_/p 2 |(M„-M)p nin - v /p 2 mn -l| 

where p mn is the cut-off ratio. This is equivalent to the 
expression that Rice, et.al. 6 , gave for cos 0'. 

iv) Other Results : One can easily show graphically 
some other known results. One is the following: if a 
mode (m,n) propagates for M=M,, then it will propagate 
for duct Mach number M 2 >M t . The proof is simple. 
From Fig. 6 the semi-minor axis of the ellipse for M 2 is 

l/p 2 , where P 2 =^1 -M \ <P 1 =^1 ~-Mf . Therefore, 

l/P 2 >l/p, = semi-minor axis of the ellipse for M,. 
Hence, a horizonal line that intersects the ellipse for 
M=M, necessarily intersects the ellipse for M=M 2 . See 
Fig. (6). 
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Another result is that as M-*l all modes will 
propagate. Again the graphical proof is very simple. 
As M-*l, P-»0 and the equation of the ellipse 
degenerates into the parabola 

1^ = 1 -2k, < 23 > 

The axis of this parabola is the kj-axis and the parabola 
intersects the kj -axis at kj = 1/2 . The parabola extends 

to infinity in the ^ direction so that any horizontal line 
intersects it. See Fig. 7. This means, by our graphical 
method of discovering propagating modes, that all 
modes propagate. 

We mention here that the ellipse of Fig. 1 remains 
unchanged for2-dimensional, 3 -dimensional rectangular, 

and annular ducts. Only the values of kj for the modes, 
which depend on the duct geometry', will change. 

Detection of Mode Radiation From a Duct 
by an External Microphone Array 

In this section we describe the mathematics of 
mode detection by use of a circular microphone array. 
The axis of this array coincides with engine axis. We 
consider a static engine here (M=0). We assume that 
we can approximate the radiation from the inlet by the 
Rayleigh formula 



where U is the axial acoustic particle velocity at the 

inlet, and R is the radiation distance from points in the 
inlet plane. We will concentrate on a frequency of 
a x BPF = a B Q , where B is the number of fan blades 
and Q is the angular velocity of the fan. Let us assume 
that we have V vanes, or a circumferentially periodic 
disturbance of period 2 tt/V radians. We know that we 
generate circumferential modes m = aB+qV where 

q = ± 1 , ±2, Let n be the radial mode number 

index of the acoustic waves that propagate and radiate 
out of the inlet. The acoustic pressure in the duct for 
frequency aBQ is 


p / = expi[aB(Qt -0)]EEA nq J m [k T (m,n)r] x 

n q 

exp{ -i[qV0 +k a (m,n)x]} (25) 


Note that m is a function of q in Eq. (25). 

The relation between u and p' for frequency aBQ 

is: 

u = expi[aB(Qt -0)] x 

aBQp 0 

EEA^k a (m,n) ^[It^m^r] x 

n q 

exp(-i[qV0 + k a (m,n)x]} (26) 

Now let us assume that the origin of the cylindrical 
polar coordinates is at x=0. Then 

u(r,0,O) = — expi[aB(Qt -0)]x 

aBQp 0 

E E A^ k a (m, n) J m [k, (m, n) r] e ’ iq ve (2 7) 

n q 

Let a microphone on a circular array of radius a from 

the x axis have cylindrical polar coordinates (a,0 / ,x). 
From a point on the inlet with coordinates (r,0), the 
radiation distance is 

R-Rq-t sin i|/ cos(0-0 / ) (28) 


where R 0 = v^a 2 +x 2 is the distance of the microphone 
from the center of the inlet and sini(/ =a/R 0 , i.e. ifi is 
the angle that the microphone makes with the engine 
axis. We are assuming that r a /R 0 «l, where r 0 is the 
inlet radius. 
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Using Eq. (27) in Eq. (24), we get 


r o 2n 


4 np / (a, 0 / ,x,t) = — expif aB[Q(t - — ) - O']) ££A nq k J (m,n)e' iqVe j /J m [k t (m,n)r] x 
P 0 c O q o O 

) i I- (aB + qV) (6 - O') +■ cos (0 - 6 ') sin i|r l rdrdG 


= — expiiaB 

Po l 


Q 


expi 


c 


- 0 / 


EE A nq k,(m ) n)e- iqve y J m [k t (m,n)r] J„ 
J n q o 


— D(0')expHaB 
Po 


Qlt- — 1 - 0 ' 


aBQr . , 
simjr 


rdr 


(29) 


where D(0') is a periodic function of 0' with first 
harmonic 2k/ V. 

From Eq. (29), we see that since 

0(0') (30) 
q 
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Fig. 1 - The ellipse of wave number vector described by Eq. (8); 

k.=k a /k,k 2 =k r /k. Corresponding to each propagating wave with a given 
kj, two axial wave numbers k 1+ and k,_ are obtained. 


Upstream 



M =0.8 


♦ Propagating 
X Decaying 


Fig. 2. - Determination of propagating and decaying waves. Each vertical line on 
the right corresponds to one circumferential mode. 
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Fig. 3. - The wave number vector at cut-off 



Fig. 4. - The diagram for finding the direction of group velocity in the k, kj -plane. 



Fig. 5. - The diagram for construction of 0', the angle of the main radiation lobe 
from an inlet, as proposed by Rice, Heidmann and Sofrin 5 . The ellipse is 
based on the duct Mach number and the flight Mach number is M^. 
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Fig. 6. - The ellipses for M=0.6 and M=0.7. Mode (m,n) propagates for M-0.6 
Therefore it propagates for M>0.6. 



Fig. 7. - The degenerate ellipse (a parabola) as M— >1. Note that as k,- -0 °, we 
have ± «. Any horizontal line intersects the parabola, i.e., all modes 
propagate. 
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